Determination of the structure of ^ Ne by full-microscopic framework 
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We perform the first quantitative analysis of the reaction cross sections of ^*~^^Ne by ^^C at 
240 MeV/nucleon, using the double-folding model (DFM) with the Melbourne g-matrix and the deformed 
projectile density calculated by the antisymmetrized molecular dynamics (AMD). To describe the tail of the last 
neutron of ^^Ne, we adopt the resonating group method (RGM) combined with AMD. The theoretical predic- 
tion excellently reproduce the measured cross sections of ^*~^^Ne with no adjustable parameters. The ground 
state properties of "^^Ne, i.e., strong deformation and a halo structure with spin-parity 3/2^, are clarified. 

PACS numbers: 21.10.Gv, 21.60.Gx, 24.10.Ht, 25.60.Dz 



Introduction. Exotic properties of nuclei in the "island of 
inversion" fascinate many experimentalists and theoreticians. 
The term "island of inversion" was first introduced by War- 
burton et al. ||Ti|] to specify the region of unstable nuclei from 
■^°Ne to '^^Mg. According to many experimental and theo- 
retical studies, it turned out that the low excitation energies 
and the large B{E2) values of the first excited states of nu- 
clei in the island indicate strong deformations i^-Q]. which 
eventually cause the melt of the neutron shell corresponding 
to the = 20 magic number (N: neutron number). To un- 
derstand these features, intruder configurations of the single- 
particle orbits have been discussed with the shell model |@], 
the generator coordinate method i^, and the Hartree-Fock 
(HF) and Hartree-Fock-Bogoliubov (HFB) methods based on 
the Skyrme or Gogny interaction lllC)l - [l3ll . In particular, the 
■^^Ne nucleus is very interesting in view of its intruder con- 
figurations and a halo structure due to strong deformations. 
Recently, a systematic investigation employing the antisym- 
metrized molecular dynamics (AMD) with the Gogny DIS 
interaction has been performed for both even and odd N nu- 
clei in the "island of inversion" ifTill . AMD was shown to give 
rather large deformations and small separation energy of the 
■^^Ne as indicated by the preceding studies mentioned above. 

Very recently, experimental studies took large steps toward 
exploring the "island of inversion", i.e., the one-neutron re- 
moval cross section (t_„ of ■^^Ne at 230 MeV/nucleon was 
measured by Nakamura et al. ifTsIl and the interaction cross 
section ctj of 28-32|^g 12q 240 MeV/nucleon was mea- 
sured by Takechi et al. |[l6l [itIi . According to the analy- 
sis of the (T_„ ifTsl [I8I - I20II . the measured large cross section 
suggested a neutron halo structure of '^^Ne with spin-parity 
3/2~. This conjecture was confirmed by the analysis of the 
Cj fl6l[l7ll . in which neutron configurations with lower partial 
waves were favored to explain the ctj of ■^^Ne. In Ref. UtIi . 
it was shown that s-wave configuration gave slightly better 
agreement with the data, which implied a possibility of even 
more drastic change of the nuclear shell structure. Note that cti 
differs from the total reaction cross section ctr by the inelas- 
tic cross section(s) due to the excitation of the projectile. For 
unstable nuclei, however, the two are almost identical, since 
very few bound excited states exist for such nuclei. Thus, cri, 
as ctr, can be assumed to represent the size of the nucleus. 



All aforementioned studies indicate the shell evolution of 
■^^Ne and its halo structure. Quantitative understanding of 
these properties of "^^Ne is, however, still under discussion. 
This is mainly due to the fact that nuclear many-body wave 
functions obtained by the high-precision structural models 
have never been directly applied to reaction calculation. As a 
first step, in Ref. IJH] we analyzed the CTj of '^^'^^Ne lfl6l[l7|l 
with the microscopic double folding model (DFM) based on 
the Melbourne 5-matrix I22ll . We adopted the mean-field 
wave functions based on a deformed Woods-Saxon potential, 
with the deformation parameter evaluated by AMD. It was 
shown that the deformation of the Ne isotopes was indeed im- 
portant to reproduce the experimental data. The agreement 
between the calculation and the data were, however, not fully 
satisfactory. The large difference in ctj between -^^Ne, '^^Ne, 
and ■^^Ne could not be explained well in particular. 



In this Letter, we directly incorporate the AMD wave func- 
tions of ^^~'^^Ne in the DFM calculation and see how the 
structural properties of these nuclei based on AMD are "ob- 
served" through (jj. For ■'^Ne, we further utilize the resonating 
group method (RGM) to give a proper behavior of the neutron 
wave function in the tail region. This is the first microscopic 
calculation of (Tj for ^^^■'■^Ne with no free parameters. The 
predicted values of ctj are validated by the comparison with 
the experimental data lfl6l [itIi . Through this study, we aim to 
determine the ground state structure of "^^Ne. 



Theoretical framework. We calculate the total reaction 
cross section cTp by DFM as in Ref. ll2lll . This model is ac- 
curate, when the projectile breakup is small. This condition is 
well satisfied for scattering analyses here, since the breakup 
cross section is quite small even for scattering of ''^Ne with 
small neutron separation energy ll20ll . A microscopic optical 
potential U between a projectile (P) and a target (T) is con- 
structed by folding the effective nucleon-nucleon (NN) in- 
teraction with the projectile and target densities, pp and prp, 
respectively. The direct (Uur) and exchange (C/ex) parts of 
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(1) 



the folding potential are obtained by 11231 l24ll 

U^RiR) = J Ppirp)PTirT)voniP^s)drpdrT., 

X WgxlPi ■s) oxp [ii^(i?)-s/A/]drp(irrp, (2) 

where s = rp — rx + -R for a position vector R of the 
center-of-mass of P from that of T. The original form of 
Ubx is a non-local function of R, but it has been localized 
in Eq. (|2|i with the local semi-classical approximation ll25ll . 
where hK{R) is the local momentum of the scattering con- 
sidered and M = ApAt^/{Ap + At) for the mass number 
Ap (At) of P (T). The validity of this localization is shown in 
Ref. ll26ll for nucleon-nucleus scattering; note that this is also 
the case with nucleus-nucleus scattering. In Eqs. ([TJ and ^ 
the direct (exchange) component of the effective NN inter- 
action, wdr (wex), is assumed to depend on the local density 
at the midpoint of the interacting nucleon pair. We adopt the 
Melbourne g-matrix as an effective NN interaction in nuclear 
medium. 

As for Pt, we use the phenomenological ^^C-density de- 
duced from the electron scattering 12711 . with unfolded the 
finite-size effect of the proton charge following the standard 
manner ll28ll . 

The projectile densities pp of ^'^^■^^Ne are calculated from 
the AMD wave functions that successfully describe the low- 
lying spectrum of Ne isotopes ifTill . The reader is directed 
to them for the details of the AMD wave function and cal- 
culation of pp. To investigate one neutron-halo nature, we 
have performed more sophisticated calculation for ■^^Ne, that 
is called AMD-RGM below. Employing the coupled-channels 
RGM type wave function. 



vI/(3iNe;3/2r) = ^ ^ {x„,Mi^™(r)«'rNe; J:)</>„} 
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(3) 



the relative wave function Xnj between the last neutron and 
the core ("^"Ne) is calculated by solving the RGM equation. 
Here the wave functions of ^°Ne are those of AMD ob- 
tained in Ref ifTill and includes the many excited states with 
positive- and negative-parity below 10 MeV in excitation en- 
ergy. Therefore, note that the weak-binding feature and the 
possible core excitation associated with the strong deforma- 
tion are properly treated in the AMD-RGM calculation. 

If one or both of the densities pp and Pt are non-spherical, 
the microscopic potential U is not spherical. It follows from 
Ref. ll2lll29ll . however, that it is sufficient to use an angular- 
averaged density in DFM, i.e., Eqs. ([T]) and (|2]i, in the present 
case. Nevertheless, deformation effects of the wave functions 
are taken into account as shown in Ref. j21]. 

Results and Discussions. We show in Fig. [T] the result of 



CTp for ^'^C, '^^He., "'"'Na, and "''Al by a ^^C target at around 
250 MeV/nucleon, compared with the experimental data Isol - 
[32I1 . For •^"Ne and ^■'Na, the original experimental data were 
measured at around 950 MeV/nucleon ll3ll[32ll but the present 



ones are corrected in 240 MeV/nucleon by Glauber calcu- 
lation 13011 . The projectile densities are phenomenological 
ones obtained by electron scattering iEtIi . The theoretical re- 
sults of (Tj^ shown in this Letter are reduced by 1.8% as in 
Ref. This fine tuning has been done to reproduce the 

mean value of the o-p of the ^^C-^^C scattering measured at 
250.8 MeV/nucleon ll30ll . Thus, the present calculation con- 
tains no free parameters except for the ^^C-^^C scattering. 
Figure [T] shows the high accuracy of the ctj^ predicted by the 
present DFM calculation for the ^°Ne, ^■^Na, and ^ ""Al projec- 
tiles. 
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FIG. 1: (Color online) Reaction cross sections for scattering of stable 
nuclei from ^'^C at 250 MeV/nucleon. The solid line represents the 
results of the present double fol ding model calculation. The experi- 



mental data are taken from Ref. I3(i - l32ll . 



TABLE I: The spin-parity and deformation parameter /? and 7 of the 
ground states of Ne isotopes calculated by AMD. 



nuclide 


28Ne 


2^Ne 


^"Ne 


3iNe 


32Ne 


r 


0+ 


1/2+ 


0+ 


3/2- 


0+ 




0.28 


0.43 


0.39 


0.41 


0.33 


7 


60° 


0° 


0° 


0° 


0° 



The structural properties of •^^^'*^Ne obtained by AMD are 
shown in Table U (spin-parity and deformation parameter /3 
and 7 defined by the Hill-Wheeler coordinate) and Fig.|2](one 
neutron separation energies Sn)- In the latter, experimental 
data issi [3 are also shown. One clearly sees that ^^^-^^Ne 
are strongly deformed and the odd-even staggering of 5„ mea- 
sured is reproduced very well. It is thus expected that the 
AMD wave functions of Ne isotopes are highly reliable. How- 
ever, for ■^^Ne, which has very small Sn (^ 250 keV), the tail 
of the wave function of the last neutron may not be described 
properly, since AMD uses a one-range Gaussian wave func- 
tion for the motion of each nucleon. 

This possible shortcomings can be overcome by using 
RGM, which generates a proper asymptotics of the last neu- 
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FIG. 2: (Color online) One neutron separation energy S„ of the Ne 
isotopes. The experimental data are taken from Ref. (33l . [34ll . 



tron. The neutron one-body density p„(r) of "^^Ne thus ob- 
tained is shown in Fig. [3] by the solid line. The results for 
■^^Ne (dashed line) without RGM are also shown for compar- 
ison. The solid line has a long tail, whereas the dashed line 
rapidly falls off at r > 6 fm. The root mean square radius 
of the density obtained by AMD-RGM (AMD) is 3.617 fm 
(3.490 fm). Although the difference between the two values 
is not so large, the density in the tail region, i.e., r > 6 fm, 
has a significant contribution to cr^ as shown below. Note 
that AMD-RGM gives 5„ = 450 keV, which is slightly larger 
than 250 keV obtained by AMD but still consistent with the 
measured value, S„ = 0.29 ± 1.64 MeV iH. 
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FIG. 3: (Color online) The neutron one-body densities of Ne. The 
solid and dashed lines represent the results of '^^Ne calculated by 
AMD-RGM and AMD, respectively. 



The AMD-RGM and AMD wave functions of ^^Ne can be 
decomposed in terms of the spherical basis components as 
shown in Table [III Compared to the AMD, the amount of 
2+ states reduces in AMD-RGM. This is due to the weak- 



TABLE II: Configurations of the ground state of '^^Ne obtained by 
AMD-RGM and AMD. 



Spherical basis 


Amplitude 








AMD-RGM 


AMD 


^°Ne(0+) 5 


5 lP3/2 


56% 


37 % 


^''Ne(2+) (> 


5 lP3/2 


24 % 


41 % 


^°Ne(2+) (> 


5 0/7/2 


9% 


12 % 


^°Ne(l-) (5 


§ lSi/2 


5 % 


5 % 


other components 


6 % 


5 % 



coupling between ■^"Ne and valence neutron. In the AMD- 
RGM model, the main component of the last neutron is 1^3/2 
coupled with the O"*" ground state of ■^"Ne. The long range 
tail of Pn(r), i.e., the halo structure, shown in Fig.[3]is due to 
this configuration of "^^Ne. Thus, reliable nuclear wave func- 
tions of Ne isotopes, including a p-wave halo nucleus ''^Ne, 
are prepared. 
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FIG. 4: (Color online) Reaction cross sections of Ne isotopes by 
^^C at 240 MeV/nucleon. The solid line represents the result of the 
present study. The result obtained by HF calculation is shown by the 
dashed line. The triangle represents the result for ^^Ne without RGM 
prescription. The experimental data are taken from Refs. 



In Fig.|4]we show by the solid line the prediction of for 
28-32|^g H- scattering at 240 MeV/nucleon. The agree- 
ment of the present result with the experimental data is excel- 
lent. If we adopt spherical HF wave functions for Ne isotopes, 
the dashed line is obtained, which significantly undershoots 
the data. More seriously, no bound-state solution is found for 
3i,32jyg The triangle shows the result for ^^Ne calculated 
with AMD wave functions. Clearly, the result cannot explain 
the very large experimental value. We conclude from these 
findings that i) ^®~'^^Ne are strongly deformed as discussed 
in Ref. ll2lll and ii) '^^Ne has a halo structure due to the last 
neutron in the orbit. Very recently, in Ref. ifssll the par- 
ing anti-halo effect was discussed for odd-even staggering of 
(Tr for ■^°~'^^Ne. It will be interesting to consider the pairing 
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effects in the present framework. To see how the strong defor- 
mations indicated by AMD affect the resuh in Ref. II35I1 will 
also be very important. 

Summary. We have performed a microscopic calculation of 
the reaction cross sections for neutron-rich Ne-isotopes sys- 
tematically. The double-folding model (DFM) with the Mel- 
bourne g-matrix and the antisymmetrized molecular dynam- 
ics (AMD) wave functions were used. AMD is a powerful 
tool that describes strongly deformed nuclei. For a loosely- 
bound nucleus ■^^Ne, the resonating group method (RGM) was 
adopted to generate a proper behavior of the wave function of 
the last neutron. The present framework reproduced the ex- 
perimental data very well with no free adjustable parameters. 
Thus, the AMD wave functions of the Ne isotopes (AMD- 
RGM for "^^Ne) were clearly validated. We concluded that 
neutron-rich Ne-isotopes are strongly deformed and ^^Ne has 
a halo structure with spin-parity 3/2~. In the near future, we 



will apply the present framework to investigate properties of 
Mg-isotopes, for further exploration of the "island of inver- 
sion". 
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